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Abstract 

Superheavy (SH) quasistable particles generated in the Early Universe could 
be responsable for Ultra High Energy Cosmic Rays (UHECR) and be a com- 
ponent of Cold Dark Matter (CDM) in the universe. These particles are likely 
to cluster in the galactic halo, so that the main part of UHECR are gamma 
rays produced in the decay of neutral pions. Charged pions are also produced 
in the same process and result in high energy electrons. We consider here the 
production of gamma rays by synchrotron emission of these electrons in the 
galactic magnetic field. The gamma ray fluxes are in the region of interest 
for some current and proposed experiments (e.g. EGRET, GLAST, MILA- 
GRO) in the energy range 0.1 — 10 4 GeV. A comparison with the existing 
upper limits at 10 5 — 10 8 GeV is also carried out. The detection of this flux of 
gamma rays would be an important signature of SH relic particles as sources 
of UHECR and would give a clue to the physics of the Early Universe. 
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I. INTRODUCTION 



The detection of cosmic rays well above 10 20 eV poses a serious challenge to our under- 
standing of the acceleration processes as well as of the propagation of these particles from 
the sources. Particles with E > 10 19 eV are thought to be of extragalactic origin since the 
galactic magnetic field is unable to confine and isotropize them. However if the sources are 
cosmological then the UHECR spectrum would have as a unique signature the presence of 
a pronounced cutoff (the so called GZK cutoff 0) at E ~ 5 • 10 19 eV, which is not observed 
(see for a recent discussion). This limits the distance to the sources of UHECR to be less 
than ~ 50 Mpc and no plausible source has been identified by current experiments within 
this distance, unless very large deflection angles are assumed, which looks incompatible with 
the current limits on the intergalactic magnetic field (see || for a review and [Q for revised 
limits in a inhomogeneous universe). 

The lack of counterparts to the observed events, together with the difficulty of the known 
acceleration mechanisms to reach the observed energies, prompted a new class of models 
(top-down models) (for a recent review see |j] and references therein) where UHECR are 
produced by the decay of grand unification massive bosons (generically named X-particles) 
with mass mx ~ 10 13 — 10 16 GeV. These X-particles can either be produced in processes 
involving topological defects (TD) or by the decay of SH quasistable (lifetime r much larger 
than the age of the Universe to) particles relics of the early universe HBHlOU in the assumption 
that they represent some fraction of the CDM. 

In the latter case, the UHECR spectrum is dominated by the halo component and hence 



it has no GZK cutoff. Signatures of this model have been discussed in [[L2|], where it was 
pointed out that if UHECR above ~ 10 19 eV are produced in the halo by the decay of 
X-particles, then the flux should be dominated by gamma rays produced by the decay of 
neutral pions from the hadronic cascade initiated by ultra high energy partons (X — > qq, 
q — > hadrons). The model is also characterized by some degree of anisotropy due to the 
asymmetric position of the sun in the galaxy [0-[13|. 
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Together with neutral pions charged pions are also generated which in turn decay into 
electrons and positrons (hereafter we shall call them both electrons) and neutrinos. 

In this paper we calculate the flux of gamma rays which is produced by synchrotron 
emission of the electrons in the galactic magnetic field when the density of SH particles in 
the halo is normalized to fit the observed UHECR flux. Our calculations show that in most 
of the cases the predicted gamma ray emission is compatible with the upper limits imposed 
by experiments like HEGRA, EASTOP, CASA-MIA and Utah-Michigan at 10 5 - 10 8 GeV 
and, in the energy range 0.1 — 10 4 GeV might be detectable by future experiments like 
GLAST. 

The paper is planned as follows: in section |H] we present a benchmark calculation where 
we estimate the magnitude of the gamma ray flux from ultra high energy electrons (UHEE) in 



the galactic magnetic field. In section III we introduce the fragmentation functions both for 



ordinary QCD and in the supersymmetric generalization proposed in fL4| . Our conclusions 
are presented in section [TV. 



II. AN ESTIMATE OF THE EFFECT 

The SH particles we are interested in are quasistable, which means that their lifetime Tx 
is much larger than the age to of the Universe. If n Q is the average density of SH particles 
with mass mj in the halo, the decay rate can be estimated as h h x° « n Q /rx- The decay of 
an X-particle typically results in the production of a quark-antiquark pair 

X -> qq 

and each quark produces a jet of hadrons (mainly Nucleons (N) and pions (ir)) with total 
energy in the jet Ej et = rjij/2. The hadron spectrum is given by the fragmentation function 
Wi(x,mx) where i = N,n and x = Ei/Ej et = lEijmx- Many discussions can be found in 
the recent literature about approximations to the fragmentation function and extrapolations 
of the low energy phenomenological expressions to the extremely large center of mass energy 
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involved in the production of UHECR by the decay of X-particles (for a recent review see 
||; a numerical approach involving the use of a QCD event generator was used in 0). We 
address this problem in the next section, while here we shall use a very simple approximation 
proposed in |L5],[16| and valid for i < 1. In this approach the spectrum of the pions, which 
represent most of the content of the jet, is given by 

= i5 x -s/»(i _ x ) 2 «i| x -8/a (i) 
dx 16 16 

The spectrum of the nucleons has the same form as in eq. (|l|) but the normalization is such 
that the ratio N/n in the jet is ~ 0.05 as suggested by the decay Z — ► hadrons at LEP. 
Moreover each pion family will take approximately 1/3 of the total pion content in the jet. 

The small N/n ratio at the production implies that on galactic distance scales the flux 
of UHECR produced by the decay of X-particles is dominated by the gamma rays produced 
in the decay of neutral pions. Therefore, in the following we shall estimate the rate n^ 1 " 
of decay, requiring that the flux of gamma rays equals the flux of UHECR at some energy. 
For the sake of simplicity we shall also assume here that the SH particles are distributed 
homogeneously in the halo whose size is Rh- The flux of gamma rays with energy E 1 at the 
Earth is then easily derived to be 

rUHE !TP \ o w l -haioRn f E ^ ,„ dN n 2 5 . halo R H ( m x \ l/2 3/2 



nW = 2x : C- dEv^P-— = -hT°— — ETjTe, (2) 

where the factor 2 in front takes into account the two jets and the factor 1/3 accounts for 
the 7r°'s only. Moreover we used dN^/dE^ = (2/mx)dN 7T /dx. 

The comparison of our prediction with the observed flux of UHECR at E = 10 10 GeV 
immediately gives 

h h /° « 1.8 ■ 10- 42 ( ™Z„ )~ 1/2 ( T^?~) cm-h-\ (3) 
x \10 14 GeVJ \100kpcJ v ; 

A fraction 2/3 of the energy in the jet is transformed into charged pions which contribute 
an electron flux through production and decay of muons. The density spectrum of these 
electrons can be calculated in the standard way (muon decay) and gives 



q e {E e )--n x K . (4) 

The ultra high energy electrons from pion decay are produced in the magnetic field of 
our galaxy, whose average strength will be denoted here by (in fiG). The synchrotron 
emission of these electrons will typically result in the production of photons with energy 

E. 2 



E - 1H0!S 'IS^J GeV. (5) 

Electrons with energy E e > 10 17 eV will then radiate photons with energy > 0.1 GeV. 
Here we estimate the flux of these photons. The synchrotron emission is the dominant 
channel of energy losses at the energies of interest, and the time for losses is so small that 
all the electron energy is very rapidly radiated in the form of gamma rays, with flux 

w - t*<*>SS - 12 ■ 1033 ( «f (dfe) v-^-vw- 

(6) 



where jj^ has been derived from eq. fl5|). The normalized flux of gamma rays can now be 
calculated using eq. (|3|) which gives the following very interesting result 

I y (E 7 ) = 2.1 ■ lO^B-^E- 7 ^ cwr 2 GeV- 1 s- 1 sr- 1 (7) 

independent on the size of the halo Rh and the mass of the X-particles mx- 

The result is also weakly dependent on the average magnetic field B^. For B^ ~ 0.1 — 1 fiG 
the flux derived from eq. (0) turns out to be detectable by EGRET in the energy range 100 
MeV< E 1 < few GeV. 

The inverse dependence of J 7 on B^ can be easily explained on the basis of the electron 
spectrum resulting from the fragmentation function: such spectrum contains most of the 
energy at high electron energies. Therefore, at fixed E 1 the larger electron energies (corre- 
spondent to smaller magnetic fields) contribute more, as far as the required electron energies 
are below mj/2. 
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The result quoted in eq. (0) changes appreciably as a function of the energy where SH 
relic particles provide the main contribution to the UHECR flux. This is more evident if we 
combine eqs. (0) and (§) to give: 

J 7 (£ 7 ) = 8.2 • 10 3 ^ HE (E UHE )Ef/ 2 HE B^% 7 /\ (8) 

where now Euhe is the energy at which we normalize the UHECR flux I^ HE (measured in 
units of cm~ 2 s~ 1 sr~ l GeV~ 1 ). From eq. (H) it results that if the normalization is carried 
out at 5 • 10 10 GeV, then the numerical factor in eq. (^) is ~ 18 times smaller. In this case, 
for ~ 1 the flux drops below the detectability level of EGRET, though it remains still in 
the detectability range of next generation satellites (e.g. GLAST). 



Let us come now to another interesting point: following the approach used in [14] and 



12], the rate of decay of SH particles in the halo may be used to normalize the same quantity 



in the extragalactic space, through the relation 

~halo ~halo 

' l X _ FCDM ^ 
Keg ^CDMpcr 

where h eg is now the rate of X particle decay in the extragalactic space, Pcd°m = 0-3 
GeVcm -3 is the energy density of CDM in the halo, Sic dm is the fraction of CDM and p cr 
is the critical density of the Universe. We shall adopt here flcDMh 2 = 0.2. If we call B^f 
(in fiG) the extragalactic magnetic field, the flux of gamma rays due to the synchrotron 
emission of UHEE can be calculated similarly to eq. (|6|) with Rh substituted by cto, where 
to = 2.06 • 10 17 /i _1 s is the age of the Universe and h is the dimensionless Hubble constant. 
Of course this treatment is valid as far as gamma rays are not absorbed on scales much 
smaller than cto, which limits this estimate to < 500 — 1000 GeV. In this assumption 
the ratio of gamma rays with energy E 1 produced in the extragalactic space relative to the 
ones produced in the halo is given by 



g(gr) = Q 23 ( Rh \ 1 (W^\ (B \ 1/4 



J 7 (£ 7 ) ' \100kpcJ V °- 2 ) \ B > 



(10) 



where we used h = 0.6. For extragalactic magnetic fields of order 10 9 — 10 10 Gauss, and a 
typical galactic field of order of pG, the two gamma ray fluxes are comparable. Unfortunately 
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only upper limits are currently available on the extragalactic magnetic field. Therefore in 
the following we shall not consider the extragalactic flux in any more detail and we shall 
concentrate our attention on the gamma ray flux from the halo. 



III. A DETAILED CALCULATION 

The simple calculations presented in the previous section suggest that the synchrotron 
emission of UHEE in the galactic magnetic field can produce an interesting flux in high 
energy gamma rays. The rate of production of electrons is normalized by requiring that the 
gamma rays generated in the decay of neutral pions saturate the flux of UHECR above some 
energy. In this section we directly use as normalization point 5 • 10 19 eV, rather than the 
10 19 eV used before, mainly inspired by two factors: 1) using more realistic fragmentation 
functions and normalizing our gamma ray fluxes at 10 19 eV it is easy to overproduce the 
UHECR at higher energy; 2) in the same conditions the flux of primary gamma rays can 
easily exceed the CASA-MIA limits at ~ 10 8 GeV. As a consequence of this more conser- 
vative approach, the fluxes that will be derived here will be appreciably smaller than the 
ones derived in the previous bentchmark calculation, and can be considered as a lower limit 
prediction. 



As first pointed out in and []T2[ the UHECR flux in this halo model is the generation 
spectrum resulting from the initial fragmentation process, hence it is dominated by gamma 
rays with respect to nucleons. 



In this section we assume, following [12], that the rate of particle production is 



parametrized as the dark matter distribution [17 



nx{R) = nr^r ( n ) 

(R/r o y[l + (R/r o rf-^ a 

where Uq is the normalization rate, R is the distance to the galactic center and r is a 
distance scale between 5 and 10 kpc. The parameters a, [3, 7 = (2,2,0) correspond to an 
isothermal profile, while a,/3, 7 = (2,3,0.2) give the best fit to observational data |I7J and 



finally a, /3, 7 = (1, 3, 1) are obtained in the numerical simulation of ref. [18 



The calculation of the spectrum of hadrons produced in a single X-particle decay is given 
by the fragmentation function for the process q —>■ hadrons. The fragmentation functions 
we use here are the MLLA limiting spectrum of ordinary QCD |L9] and its supersymmetric 
generalization given in []nj (hereafter SUSY-QCD). In the first case the fraction of pions 
produced in a single event is taken as ~ 1, while in the second case we take, according 
with ||14|| , /V ~ 0.5. The pion fragmentation functions, indicated by W n (x,m x ) are then 



normalized as 

/ dxxW n (x,m x ) = f n (12) 
Jo 

where x = 2E 7T /m x . Approximately 1/3 of is carried away by neutral pions while a 
fraction 2/3 goes into charged pions. Each X-particle will produce two jets with energy 
mx/2 each. 

The choice of the fragmentation function and its normalization are delicate points and 
deserve some additional comments: both the MLLA fragmentation functions considered here 
(QCD and SUSY-QCD) are still valid for x -C 1, though they represent an improvement 
with respect to the simple approximation in eq. (|I|). This implies a lack of consistency in 
the way we normalize the fragmentation function through eq. (|T^) since the integral extends 
to a range of values of x where the fragmantation function is not accurate. Recently this 
problem has been treated more carefully in |7| and several new interesting results have been 
obtained. In [0] the HERWIG QCD event generator was used to investigate the structure 
of the hadronic cascade at large values of rax- in this approach the normalization of the 
fragmentation function can be carried out self consistently. More interestingly however, the 
authors of also find appreciable differences between the fragmentation functions produced 
by their montecarlo tecnique and the functions generally used in these calculations as well 
as in the present work. The differences, as could be expected, vanish for mx < 10 3 GeV, 
where the MLLA approximation is known to work well. A striking difference between the 
montecarlo approach and the MLLA approximation is that the fragmentation functions for 
protons and pions (and therefore neutrinos and gamma rays) have different shapes for large 
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values of x (while we assumed here that the shape of the two functions is the same but 
with a different normalization). Moreover, for 0.2 < x < 0.4 the authors of J7| seem to 
obtain equal abundancies of nucleons and gamma rays, which is a very different result from 
the usual assumption that the gamma rays are a factor ~ 10 more abundant than nucleons 
in the hadronic cascade. The combination of different normalization and different shape 
of the fragmentation functions suggested in can affect the conclusions presented in this 
paper and more in general the calculations of particle fluxes in top-down models of UHECR. 
Unfortunately, since the results in |7fl are purely numerical it is not easy to envision in 
what direction the results presented here would change, therefore, we assume here that the 
fragmentation process is well represented in the MLLA approach introduced above, stressing 
that a detailed study of the problem would certainly be useful. 

Since the electrons radiate through synchrotron emission, a model for the galactic mag- 
netic field is also needed. We parametrize the magnetic field in the form 

B(r,z) = B (r)exp(-z/z ) (13) 

where B (r) = 3(r Q /r)//G when the radial distance r in the disc is larger than 4 kpc and 
B (r) =const for r < 4 kpc (e.g. |2"Dfl ). The parameter z is fixed at 0.5 kpc and r = 8.5 kpc 
is the distance of the sun relative to the galactic center, z is the height above the galactic 
disc. 

According to [12[] the flux of gamma rays with energy produced in the ir° decay is 



given by 

= £ dcos6\cos6\ f™ x(e) drn h x (R) (14) 



where r max {6) is easily calculated to be 

r m ax(0) = r G cos6 + \J R 2 h - r%sin 2 9 (15) 

and Rh is the size of the halo. The angle 9 is measured with respect to the direction of the 
galactic center. The gamma ray spectrum at the production is given by [12[ 



4 /"i dx 

W ( E ^) = / —WXx,m x ) (16) 

17lx J2E-y/m x X 

with W 7V o(x,mx) = (l/3)W / 7r (x, m x )- 

The comparison of the observed UHECR flux at E = 5 • 10 10 GeV with the flux of 
gamma rays from eqs. (|I^|T6l) at the same energy gives the normalizing value of for each 
specific choice of the values of the parameters. We plot our fluxes of UHE gamma rays at 
the normalization point in Fig. 1. 

The spectrum of UHEE is the convolution of the pion and muon spectra and is given by 



W e {E e ) = l2El-P^- I ^ <l4 ^W^famx), (17) 

m% - m^mx J E e E* E„ J2E u / mx x 



r2r 2 E^/mx fa 



Ez 

where r = m^/m^. The spectrum of electrons produced by a single muon is peaked at 
E e fa (2/3)E tl . In the simple assumption that all electrons are produced with this peak 
energy the above electron spectrum becomes 

Q T7? 2 r3r 2 E e /m x At 

W e (E e ) = / —W^(x,m x ). (18) 

m x - J3E e /m x x 

As pointed out in the previous section the UHEE radiate by synchrotron emission in the 
galactic magnetic field in the gamma ray energy range, according with eq. fl5|). The emis- 
sivity per unit volume in the form of gamma rays with energy E 1 at the position at distance 
R from the galactic center can be written as 

E dE 

q,(E,,r,z)=n x (R)W e (E e )-^—^ (19) 

where E e and E 1 are related through eq. (|5|) but the magnetic field is now a function of 
the position in the galaxy. Hence the same E 1 will correspond to different E e in different 
positions in the galaxy. 

Since the magnetic field is symmetric with respect to the disc while the dark matter profile 
is spherically symmetric, the integration over lines of sight requires the explicit integration 
over both 9 and <f), so that the flux of gamma rays at the Earth, per unit time, surface, 
energy and solid angle is 
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1 /-27T rl rr m ax{d) 

Ij(E 7 ) = j d(j) j dcosO \cos6\ J drq y (E yj r J z), (20) 

where r max {6) is again given by eq. (]T5|). 

The total observed gamma ray flux at energy E 7 , J* oi (i? 7 ), is the sum of the synchrotron 
gamma ray flux [eq. (ED])] and the primary gamma ray flux [eq. flU)]. 

The results of the calculations are plotted in Figs. 2 and 3, where we used tq = 10 
kpc and z = 0.5 kpc. Although the calculation has been carried out for the three sets 
of values of the a, (3, 7 parameters, the results are quite insensitive to the specific model 
of the halo distribution of dark matter, and the related curves would be indistiguishable. 
Therefore we plotted there only the case (a,/?, 7) = (2,2,0) for reference. In Fig. 2 we 
plotted El f - 75 I t ° t {E 1 ) in order to amplify the spectral differences. The two thick lines refer 
to nix — 10 14 GeV, while the two thin lines are for nix — 10 13 GeV. The solid curves are 
obtained using the SUSY-QCD fragmentation function while the dashed lines are calculated 
with the ordinary QCD fragmentation function. For comparison we also plotted by a dash- 



dotted line the extragalactic diffuse gamma ray background [28]. All our fluxes are below 
this limit in the energy range accessible to current experiments. The points with arrows 
are the upper limits on the gamma ray flux from the HEGRA [21]], Utah-Michigan p2 



EAS-TOP H| and CASA-MIA ||] experiments. 

Clearly the stronger constraints on the gamma ray flux from the halo are imposed by the 
CASA-MIA and EAS-TOP experiments. In most of the cases considered in this paper the 
fluxes are below these limits, with the exception of the case nix = 10 13 GeV with SUSY- 
QCD fragmentation function (thin solid line). Note that in this case the overproduction of 
gamma rays at ~ 10 7 — 10 s GeV is mainly due to the primary gamma rays from 7r° decay, 
rather than to the synchrotron flux. 

In Fig. 3 the gamma ray fluxes have been blown up in the energy range between 100 
MeV and 10 4 GeV (lines labelled as before) and compared with the sensitivity of some of 
the present or planned gamma ray experiments in the same energy range (e.g. EGRET, 
GLAST, MILAGRO). 
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It could be particularly interesting to look at the gamma ray signal in the region above 
~ 500 GeV, where the contamination due to the isotropic diffuse gamma ray background, 
believed to be of extragalactic origin, should be reduced due to absorption of gamma rays 
on the infrared background. Our calculations show an integrated flux above 500 GeV at the 
level of ~ 10~ 12 photons cm _2 s _1 sr _1 in our best case scenario. 

There is clearly a dependence of the results on mx, introduced by the fragmentation 
functions, more complicated than the simple power law adopted in section |T|. The difference 
is particularly evident on the tail of the spectra where the highest energy electrons (with 
energy close to ~ mx/2) are probed. 

We also checked that our predictions do not depend appreciably on the typical scale 
height z of the magnetic field (eq. fll3|)). 

A few comments about the results presented in the figures are needed: the factors 4-7T in 
eqs. ([14]) and ( p0[) scale the fluxes to the unit solid angle; this is the correct procedure only for 
an isotropic distribution of the arrival directions of the gamma rays. Due to the anisotropic 
distribution of the dark matter halo and of the magnetic field of the Galaxy as seen from the 
sun, neither the synchrotron flux nor the primary gamma ray flux are isotropic. In fact our 
preliminary results of the anisotropy for the synchrotron gamma rays |27J shows a pattern 
of arrival directions which is quite different from the one of primary gamma rays, calculated 
in 



12] and [11]. In this sense the fluxes per unit solid angle plotted in the figures must 



be intended as indicative of the 'real' fluxes, but more detailed calculations are definitely 
needed in this direction. For the same reason, summing up the contribution of primary and 
synchrotron gamma rays as done above is not self consistent for fluxes anisotropic in different 
ways. A peculiar finding in the anisotropy pattern could be important to discriminate the 
synchrotron gamma rays produced by SH particles from other background components and 
could improve the possibility to detect the signal discussed here. 
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IV. CONCLUSIONS 



We calculated the flux of gamma radiation produced through synchrotron emission in 
the galactic magnetic field from ultra-high energy electrons generated as decay product of 
charged pions in the hadronic cascade coming from the decay of X-particles, when they 
cluster in the galactic halo in the form of SH (mx ~ 10 13 — 10 14 GeV) quasistable particles. 

If the magnetic field in the galactic halo is Bh a io ~ 0.1 — lfiG, UHEE with energy between 
10 17 eV and ~ mj/2 very rapidly loose all their energy in the form of gamma rays mostly 
with energies E 7 < 10 6 — 10 8 GeV. The resulting diffuse gamma ray emission above 100 
MeV is of the order of 10~ 8 photons cm~ 2 s^ 1 sr~ l in the best case scenario. This result 
is not strongly dependent on the mass mx of the SH particles (though some dependence 
is introduced due to the fragmentation functions used here) nor on the geometry of the 
galactic dark matter halo, but depends critically on the energy where the contribution of SH 
relic particles to UHECR becomes dominant (section p|) . The shape of the fragmentation 
functions at very high energy can affect the result in several ways, and an accurate procedure 
of extrapolation of the fragmentation functions tested in accelerator experiments up to 
extremely large values of mx is needed, as stressed in [|7| . In fact in |7[ appreciable differences 
were found respect to the MLLA approximation used here, when the fragmentation functions 
are determined by using the HERWIG QCD event generator. 

In section H we also estimated the extragalactic flux of gamma rays due to the same 
physical process, and found that, for typical extragalactic magnetic fields of ~ 10~ 10 — 10 -9 
Gauss the extragalactic and the halo fluxes are comparable, if SH particles provide a fraction 
QcDAdh 2 = 0.2 of the CDM and the SH particles in the halo explain the UEHCR flux (the 
galactic and extragalactic densities of SH particles are related through eq. @). 

It is interesting to note that the detection of a residual gamma ray emission from the 
halo, after the subtraction of the contribution of pp collisions, Inverse Compton scattering, 
electron bremsstrahlung and of the extragalactic diffuse gamma ray background was recently 
claimed. This subtraction procedure is quite difficult and very model dependent; however 
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the analysis in |25| and [26| show that the residual emission is at the level of ~ 10 7 photons 



cm~ 2 s~ 1 sr~ l above 1 GeV with a spectrum E~ a and a ~ 1.8. This flux, if confirmed, is 
larger than what is found in the present calculation with the normalization at the UHECR 
flux carried out at 5 ■ 10 10 GeV, but the power index a ~ 1.8 is close to 1.75 in eq. (0). In 
the more realistic calculations illustrated in section [TT] the spectrum cannot be approximated 
by a unique power law but in the energy region between 100 MeV and a few GeV the slope 
of the spectra in the log-log plot is around 1.8 — 1.9. 

The origin of the residual emission is still under debate and the possibility of an under- 
estimate of the contribution of inverse Compton scattering is an open possibility, so that 
the actual residual emission (if any) could be smaller than what previously claimed. In this 
prospective it is interesting to study the contribution of SH particles in the halo to this 
emission as an open intriguing possibility requiring further investigation, though probably 
only next generation gamma ray detectors can shed a new light on the problem. 

A particularly interesting aspect is the anisotropy of the gamma ray emission expected 



in this model |27| . The anisotropy is the result of both the asymmetric configuration of the 
magnetic field in the galaxy as seen from the Earth, and the asymmetric position of the 
solar system in the galaxy (this also produces some anisotropy in the arrival directions of 
UHECR (HHI). 



The upper limits reported in Fig. 2 are due to the HEGRA Utah-Michigan [22 



EASTOP [^] and CASA-MIA J24| experiments. The stronger constraints to the model of SH 



particles in the halo come from the CASA-MIA and EASTOP experiments: the calculated 
fluxes are compatible with the experimental limits with the exception of the case obtained 
for nix = 10 13 GeV and the SUSY-QCD fragmentation function. In this case the exceeding 
limit is mainly due to the primary gamma rays from 7r° decay rather that to the synchrotron 
gamma rays. 

We propose that it could be of particular interest to study the gamma ray emission 
from the halo in the energy range E 1 > 500 — 1000 GeV, for two reasons: a) the diffuse 
extragalactic gamma ray background is expected to be less pronounced in this energy range 
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because of absorption of gamma rays in the intergalactic medium on the infrared photon 
background; b) the spectrum of the gamma radiation from SH relic particles is predicted to 
have a spectrum E~ a with a ~ 1.7 — 2, substantially flatter than the gamma ray emission 
from pp collisions (expected to have the same spectrum as the cosmic rays) which should 
dominate over other mechanisms at these high energies. 
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Figure Captions 



Fig. 1: Flux of UHE gamma rays from 7r° decay for the case of SH particles in the halo. 
The thick lines are for m>x = 10 14 GeV and the thin lines are for mx = 10 13 GeV. The 



solid lines are obtained with the SUSY-QCD fragmentation function |14|, while the dashed 
lines are for the ordinary QCD fragmentation function ||19|| . The fluxes are multiplied by 
E 3 . Data points are from ref. [[!]. 

Fig. 2: Total flux of gamma rays as sum of the synchrotron emission of UHEE in the 
galactic magnetic field and the contribution of the decay of neutral pions. Lines are labelled 
as in Fig. 1. Fluxes are multiplied by E^ 75 to emphasize the spectral differences among 
the curves. Upper limits are from the HEGRA 0, Utah- Michigan @, EASTOP @ and 



CASA-MIA [^3] experiments. The dash dotted line is the extragalactic diffuse gamma ray 
background from EGRET 

Fig. 3: Flux of gamma rays (multiplied by E^) for the same cases as in Fig. 2 but 
limited to the energy range 100 MeV < E 1 < 10 4 GeV. Lines are labelled as in Fig. 1. The 
sensitivity of the EGRET, GLAST and MILAGRO experiments are drawn. 
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